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ABSTRACT 

We have used far-infrared data from IRAS, ISO, SWIRE, SCUBA and MAMBO 
to constrain statistically the mean far-infrared luminosities of quasars. Our quasar 
compilation at redshifts < z < 6.5 and /-band luminosities —20 < Jab < —32 is the 
first to distinguish evolution from quasar luminosity dependence in such a study. We 
carefully cross-calibrate IRAS against Spitzer and ISO, finding evidence that IRAS 
100 fim fluxes at < 1 Jy are overestimated by ~ 30%. We find evidence for a correlation 
between star formation in quasar hosts and the quasar optical luminosities, varying as 
SFR oc LQ p 4 t 4±0 07 at any fixed redshift below z — 2. We also find evidence for evolution 
of the mean star formation rate in quasar host galaxies, scaling as (l + z) 16 * 3 at z < 2 
for any fixed quasar /-band absolute magnitude fainter than —28. We find no evidence 
for any correlation between star formation rate and black hole mass at 0.5 < z < 4. 
Our data are consistent with feedback from black hole accretion regulating stellar 
mass assembly at all redshifts. 

Key words: submillimetre - infrared: galaxies - galaxies: high-redshift - galaxies: 
active - quasars: general. 



1 INTRODUCTION 

The interaction between the accretion process into a super- 
massive black hole residing at the centre of an active nu- 
cleus and star formation in the host galaxy is fundamen- 
tal in regulating both galaxy evolution and the growth of 
the black hole. In order to understand the link between the 
, two processes and assess the possibility of the two occur- 
ing concomitantly, it is important to quantify and constrain 
the star formation activity in quasar host galaxies. This is, 
however, a difficult task, as the star formation signature can 
be diluted by the strong AGN emission, especially at short 
(e.g. UV/optical) wavelengths. Emission emerging from star 
formation activity should, therefore, be looked for in the far- 
infrared (FIR), where the contribution of the AGN (in the 
form of thermal emission from dust) should be less impor- 
tant. 

Combined AGN studies with IRAS and ISO already es- 
tablised the presence of strong FIR emission in quasars (see 
e.g. Verma et al. 2005 and references therein) . According to 
some models, this radiation might be explained as the emis- 
sion of dust distributed in a "cloud" around the central en- 
gine with a 0.5 kpc radius (Siebenmorgen et al. 2004). How- 
ever, other models estimate dusty tori extending to several 
kpc (e.g. Fritz et al. 2006, Hatziminaoglou et al. 2008) with 
additional very large covering factors (> 90%). Spitzer IRS 
spectroscopy, revealed in addition to the FIR emission, the 
presence of PAH features in the mid-IR spectra of optically- 



selected quasars (Schweitzer et al. 2006), that are difficult 
to reproduce in models assuming dust heated by the hard 
AGN photons. A number of other arguments including the 
likely evaporation of PAH features in the presence of AGN 
emission and in the absence of high column densities (a nec- 
essary condition for the AGN to be able to heat the dust 
at such large distances) or the simultaneous presence of star 
formation evidence at other wavelengths suggest that the 
FIR in quasar host galaxies is more likely to be a tracer of 
star formation. 

Even though Spitzer observations have increased the 
number of FIR detections in most low redshift quasars (e.g. 
Schweitzer et al. 2006) and their analysis consistently points 
toward star formation driven FIR emission, the number of 
FIR detected quasars is still low. In preparation for the Her- 
schel ATLAS (Astrophysical Terahertz Large Area Survey), 
an Open Time ~ 500 deg 2 blank-field survey, and in prepara- 
tion for targeted Herschel surveys of AGN, we need the best 
possible estimates for the quasar fluxes in Herschel bands. 
Most Sloan Digitized Sky Survey (SDSS) quasars are not 
detected individually in the Spitzer SWIRE Legacy Sur- 
vey 70 /im and 160 /im data (Lonsdale et al. 2003, 2004) 
in the SWIRE-SDSS overlap region (e.g. Hatziminaoglou et 
al. 2005, 2008). Submm and mm- wave observations of z ~ 2 
and z > 4 quasars have yielded only a small number of di- 
rect detections (e.g. Omont et al. 1996, Omont et al. 2001, 
Carilli et al. 2001, Isaak et al. 2002, Priddey et al. 2003a, 
2003b, Omont et al. 2003, Robson et al. 2004, Beelen et al. 
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2006, Petric et al. 2006, Wang et al. 2007). IRAS and ISO 
detected just over half of the Palomar- Green quasar sample 
at 60 /im (e.g. Sanders et al. 1989, Haas et al. 2000, 2003). 

In this paper, we will use stacking analyses to constrain 
the mean far-infrared luminosities of quasars, selected over a 
very wide range in redshift and absolute /-band magnitude. 
Our quasar compilation spans enough of the J-magnitude- 
redshift plane to be able to distinguish evolution from quasar 
luminosity dependence, which would be impossible in a sin- 
gle /-magnitude- limited quasar sample. Previous authors 
have stacked quasar fluxes at one (usually submm) wave- 
length, but we will combine a large body of multi-wavelength 
far-infrared, submm and mm-wave quasar photometry using 
an assumed common SED. We will then show that our con- 
clusions are robust to reasonable choices of SED. We follow 
SDSS in adopting the concordance cosmology of Hm = 0.3, 
Qa = 0.7, Hq — 70 km/s/Mpc, and in assuming an opti- 
cal spectral index of din S v /dlnv = —0.5 for the quasars. 
In the far-infrared we assume an M82 SED shape from Efs- 
tathiou, Rowan-Robinson & Siebenmorgen 2000, unless oth- 
erwise stated. 



2 METHODOLOGY 

2.1 Sample selection 

Figure \T\ shows the absolute I magnitudes of the SDSS 
quasars with SWIRE 70 /im and/or 160 /im coverage, against 
redshift. The SWIRE data was retrieved on 17th August 
2007 and comprises version 2 products in the ELAIS Nl 
and ELAIS N2 fields, and version 3 products in the Lock- 
man Hole field. There are 281 DR5 SDSS quasars (Adelman- 
McCarthy et al. 2007) in the SWIRE fields with 70 /zm 
and/or 160 /im coverage. Of these, 264 have 70 /im cover- 
age and 261 have 160 (im coverage. ELAIS Nl is only partly 
covered by SDSS. There is no SDSS data for the Southern 
SWIRE fields XMM-LSS, ELAIS SI or CDF-S, though this 
data was retrieved for calibration (see below); the version 
numbers are 2, 3 and 3 respectively. Figure [T] also shows 
for comparison the Palomar-Green sample with far-infrared 
photometry from IRAS and ISO (Sanders et al. 1989, Haas 
et al. 2000, 2003), in which non-detections have been re- 
measured using the SCANPI IRAS fluxes discussed below. 
Figure [1] also shows a compilation of the quasars observed 
at 850 /im and 1200 pm (Omont et al. 1996, Omont et al. 
2001, Carilli et al. 2001, Isaak et al. 2002, Omont et al. 
2003, Priddey et al. 2003a, 2003b, Robson et al. 2004, Wang 
et al. 2007). The addition of the SDSS quasars to these data 
sets greatly improves the coverage of the optical luminosity- 
redshift plane. This will make it possible to make the first 
constraints on the evolution and luminosity dependence of 
star formation in quasar hosts. The right-hand panel of fig- 
ure [1] demonstrates the far-infrared luminosities probed by 
the various multi-wavelength data sets. 

2.2 SWIRE photometry 

The SWIRE 70 /im and 160 /im images are supplied cali- 
brated to surface brightness units of MJy/sr. We intend to 
measure the total mean point source flux of our targets, 
rather than the mean surface brightness at their locations, 



Predicted conversion Emprical conversion 
70 /im 11.44 14.76 
160 /nn 59.76 81.14 

Table 1. Point source fluxes in mjy for a source in the SWIRE 
maps with a central surface brightness of 1 MJy/sr. 

so we need to convert the units. We first note that a Gaus- 
sian point source with a peak flux of unity will have a total 
flux of F = 7r#jl WHM /(4 In 2), where #fwhm is the full- width 
half maximum of the point spread function. If #fwhm is 
measured in arcseconds, then this also gives the conversion 
between total point source flux and the flux per square arc- 
second at the peak. 

Using (9fwhm = 1.22A/Z), where A is the observed wave- 
length and D — 85 cm is the diameter of the telescope's 
primary mirror, we predicted the conversion between point 
source flux and surface brightness. This is shown in table [T] 
We also derived an empirical conversion obtained by com- 
paring the SWIRE point source catalogue fluxes with the 
background-subtracted map fluxes measured at the posi- 
tions of the SWIRE sources. This comparison is shown in 
figure [5] The empirical conversion is some 29-36% higher 
than the theoretical prediction, which may be due to the 
finite size of the map pixels, or a non-Gaussian shape to the 
point spread function shape (e.g. more resembling an Airy 
function) . We adopt the empirical conversion in the analysis 
below. 

Although the maps are supplied with foreground DC 
offsets removed, we opted to subtract the mean flux level 
from each image. This ensures that the average total flux 
contribution is identically zero from each point source not 
associated with the SDSS quasars. The stacking methodol- 
ogy is to measure the point source fluxes from the SWIRE 
images at the positions of the SDSS quasars, and search 
for a significant deviation from zero flux. Photometric er- 
rors were estimated by measuring the standard deviation of 
the SWIRE images. The SDSS quasar astrometric errors are 
negligible compared to the SWIRE 70-160 /xm beams. 

2.3 IRAS photometry 

SCANPI IRAS flux estimates have recently been adapted 
to report negative AMP flux estimates, where acceptable 
fits are available. This would make this estimator suitable 
to stacking analyses, but we found positive fluxes on aver- 
age reported at positions randomly offset from our targets 
(as noted also by e.g. Morel et al. 2001). This was found to 
be due to the SCANPI algorithm allowing the source posi- 
tion to vary, so that the fits tended to gravitate to nearby 
noise features when the signal-to-noise is very low. There- 
fore, we re-fit the coadded SCANPI profiles allowing only the 
point source flux to vary, using the appropriate point source 
response function for the predominant scan direction (A. 
Alexov, priv. comm.). We used the background-subtracted 
median coadded scans. At 60 /im, we subtracted an addi- 
tional background estimate obtained by fitting a Gaussian to 
the off-source data histogram for each coadded scan (where 
"off-source" is defined as where the template is < 1% of its 
maximum value) . At 100 /im, the stronger baseline drifts vis- 
ible in the coadded scans suggested a more local background 
subtraction. We defined the 100 pm on-source width to be 
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Figure 1. (Left) Absolute I magnitudes of quasars in this paper, as a function of redshift. Filled circles represent SDSS quasars with 
SWIRE coverage; open circles represent Palomar-Green quasars with IRAS and B-band data; diamonds represent quasars observed 
at 850 /im; open squares represent quasars observed at 1200 fim (see text for details). Note that adding the SDSS quasars to the other 
samples greatly improves the coverage of the optical luminosity - redshift plane. (Right) Predictions for the far-infrared and submm fluxes 
of an M82 SED normalised to vL v = 1O 12 L0 at 100 /im. Note that the stacked signal from SDSS quasars at z ~ 2 are well-matched in 
far-infrared luminosity to the expected stacked signal from Palomar-Green quasars at z ~ 0.5. 
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Figure 2. Conversion between map surface brightness units and point source total fluxes, using the SWIRE catalogues (not restricted 
to SWIRE quasars). The 70 /im data are shown to the left, and the 160 /an data to the right. The dashed (red) line shows the predicted 
conversion discussed in the text, and the full (blue) line gives the empirical conversion. Numerical values for the conversions are given in 
tabled] Symbols: + ELAIS Nl, * ELAIS N2, o CDFS, A ELAIS SI, □ Lockman Hole, x XMM-LSS. 



where the template is > 0.1% of its maximum, and sub- 
tracted the mean of the coadded scan data in one template- 
width either side of the target. Note that our flux calibra- 
tions discussed below were found to depend weakly on the 
sky subtraction algorithm, but were not independent of it. 

We tested the flux calibration of stacked fluxes us- 
ing 70 /im sources selected from the Spitzer legacy surveys 
SWIRE and Formation and Evolution of Planetary Systems 
(FEPS, e.g. Meyer et al. 2004, Hillenbrand et al. 2008), and 
using 90 /im sources selected from the European Large Area 
ISO Survey (ELAIS, Rowan-Robinson et al. 2004). Note that 



the SWIRE and ELAIS surveys were conducted in regions 
of low cirrus; the FEPS survey, while having fewer sources, 
is more widely-distributed. We selected all SWIRE 70 /im 
sources in the flux range 12-300mJy, then starting from the 
brightest, we rejected any source closer than 30' to any se- 
lected source. This ensured that both our source fluxes and 
their background estimates are statistically independent. We 
adopted the same selection procedure for ELAIS. For FEPS, 
we selected all sources with 70 /im detections at 5a or above. 
For FEPS, SWIRE and ELAIS sources, we extracted IRAS 
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Figure 3. The left panel shows the unweighted mean IRAS 60 (im AMP SCANPI fluxes for Spitzer 70 /an SWIRE sources (circles) 
and FEPS sources (squares). The errors are the usual ojyN estimate of the noise on the mean, except for the zero-flux bins where 
the signal and noise are estimated from a Gaussian fit to the measurements as discussed in the text. The full line shows the 1:1 
relation, and the dashed lines show the expected range of colours at 0.5 < z < 1.5 for M82, Arp 220 and a cirrus spectrum, from 
Efstathiou et al. 2000. Colours were calculated using the respective wavelength-dependent system responses, following the prescription 
in http://spider.ipac.caltech.edu/staff/lord/MIPS/MIPS.html The redshift dependences of this colour for these SEDs are shown in 
the right panel, with Arp 220 as short dashes, M82 as long dashes, and cirrus as dash-dot. 




Figure 4. The left panel shows the unweighted mean IRAS 100 fim stacked fluxes, as a function of ELAIS ISO 90 fim fluxes. Errors are 
as in figure [3] The dashed lines show the mean IRAS-ISO calibration offset found by Heraudeau et al. (2002), and the RMS of their 
sample. The right panel shows the expected IRAS-ISO colours expected for the SEDs in figure [3] 



fluxes in a 3 x 3 grid centred on the target, with grid posi- 
tions separated by 20'. 

To estimate the noise in each SCANPI flux estimate, 
we tried two approaches: firstly, estimating the variance in 
the coadded scan from Gaussian fits to the data histogram, 
then propagating the noise in the best-fit point source am- 
plitude assuming the data points are statistically indepen- 
dent; secondly, measuring the variance of the point source 
amplitude estimates in the eight offset sky positions. The 
latter should give signal-to-noise histograms with zero mean 
and unit variance, while for the former the histogram of 



(Siras — SspitzcO/JViRAS (where S and are signal and noise 
respectively) should also have zero mean and unit variance. 
The SCANPI data points are not statistically independent, 
and we found that propagating the noise led to error esti- 
mates too small by a factor of around 3.6 at 60 fim, and 5.7 at 
100 fim. The offset positions had consistent noise estimates, 
but have the disadvantage that the estimates are not local 
to the target. A visual inspection of the IRAS maps around 
our Palomar-Green targets (A. Alexov, priv. comm.) sug- 
gested problematic cirrus structure near several targets. We 
therefore adopted the propagated noise estimate, scaled by 
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a factor of 3.7 (5.7) at 60 /mi (100 (im) to account for corre- 
lations between the SCANPI data points. Cirrus structure 
can in principle be alleviated with a matched filter tuned 
to the power spectrum of the background (e.g. Vio et al. 
2002), though in this case the IRAS detector responsivity 
variations would present a significant complication. 

We found the SCANPI 1002 coadded scans (median 
combination of IRAS scans) to give the best signal-to- 
noise. Figures [3] and [4] show the unweighted mean average 
IRAS fluxes of the fainter Spitzer and ISO sources in broad 
flux bins. This stacking methodology appears to give unbi- 
ased estimates even at the faintest 70 /jm fluxes tested, e.g. 
~ 40 mjy at 70 /mi. This is much fainter than the mean 
60 /im fluxes of Palomar-Green quasars, which we estimate 
to be 201±35 mjy. For the offset positions, we estimated the 
signal and noise from a Gaussian fit to the histogram of the 
measurements, in order to avoid serendipitous sources. The 
stacked 60 /im fluxes at random offset positions were con- 
sistent with zero (e.g. 3.0 ± 1.5 mjy for SWIRE, see figure 

At 100 /jm, the stacked flux at blank-field positions 
offset from SWIRE sources was —6.7 ± 6.3 mjy. However, 
we found evidence for flux calibration discrepancy between 
ELAIS and our 100 /im stacked fluxes (figure [4}. Heraudeau 
et al. 2002 found that the ELAIS 90 /im point source fluxes 
were lower than IRAS FSC 100 /im fluxes by a factor 
0.76 ±0.17, after taking into account the colour corrections. 
They attributed this to systematic IRAS overestimates in 
the faintest catalogues fluxes, citing Moshir et al. 1992; how- 
ever, we find no evidence for such a calibration offset in our 
60 /jm IRAS stacks. We compared the ELAIS and SWIRE 
flux calibrations by interpolating the 70 /im and 160 /mi 
fluxes to obtain 90 /mi flux estimates, and found a 1:1 cor- 
relation with the ELAIS fluxes, so we think it unlikely the 
ELAIS flux calibration is at fault. We also calculated the 
100/mi:90/mi colour for our three SED models (figure [4}, 
and found this too small an effect to account for the dis- 
crepancy. A similar discrepancy has been noted by Jeong 
et al. (2007), comparing IRAS 100 /im and AKARI 90 /mi 
fluxes. We have therefore chosen to adopt the Heraudeau 
et al. correction factor of 0.76 to our stacked 100 pm IRAS 
fluxes. In the appendix, we present our SCANPI photometry 
of PC quasars and discuss the problematic cases. 



3 RESULTS 

3.1 Stacked flux results 

Figure [5] shows the SWIRE 70-160 /im fluxes for the SDSS 
quasars as a function of redshift and absolute magnitude. Er- 
ror bars have been suppressed for clarity, except for a single 
representative example in each panel. The SDSS quasar pop- 
ulation appears bimodal or at least with a significant skew 
to bright fluxes for a minority of objects, with a small num- 
ber of far-IR-loud objects at both 70 (im and 160 /mi. For the 
bulk of the population, there also appears to be a significant 
positive flux at the positions of the SDSS quasars at redshifts 
z < 3, and at most absolute magnitudes. Of the quasars at 
0.5 < z < 1.5, 6/114 have 70 /mi fluxes above 20mJy, and 
6/113 have 160 /mi fluxes above 70mJy. In the 1.5 < z < 2.5 
redshift bin, 3/86 have 70 /mi fluxes above 20mJy, and 3/84 



have 160 /im fluxes above 70 mjy. The green curves show 
the predicted flux for an M82 SED, normalised to 10 11 Lq, 
1O 12 L and 10 13 L Q . 

To test whether the skewed distribution is caused by a 
distinct population of far-IR-loud quasars, or whether un- 
related far-IR-bright foreground companion galaxies are re- 
sponsible, we compared the fluxes at the positions of SDSS 
quasars with those in the map as a whole, following the 
methodology of Serjeant et al. (2004). This test is prefer- 
able to Monte Carlo randomizations of the quasar positions, 
since it compares the quasar fluxes with the entire distribu- 
tion of map fluxes, rather than a randomly-selected subset. 
The results are shown in figure H3 The observed skewness in 
the quasar far-infrared fluxes is therefore not due to a patho- 
logical distribution of map fluxes. In table [2] we present the 
average far-infrared fluxes for the SDSS quasars in redshift 
bins, both with and without the high-flux population. 

It is clear from the curves in figure [5] that although the 
SDSS quasars span a narrow range of far-infrared fluxes, 
they span a wide range of far-infrared luminosities. The 
right-hand panel in figure [T] demonstrates that the stacked 
160 /im signal from z ~ 2 quasars in SDSS covers compa- 
rable far-infrared luminosities to the expected stacked sig- 
nal (a few times fainter than the IRAS PSC) from z ~ 
0.5 Palomar-Green quasars, for the assumed M82 SED. 
Also, it is clear that combining the Palomar-Green, SDSS 
and SCUBA quasars gives enough coverage of the optical 
luminosity-redshift plane to remove the degeneracy between 
trends in luminosity and in redshift (essentially Malmquist 
bias in its original sense, Malmquist 1924; see also Teerikorpi 
1984). 



3.2 Stacked luminosity results 

It is not obvious what the best stacking statistic is where a 
large number of the sample have high signal-to-noise direct 
detections, as is the case with our quasars. Variance weight- 
ing the stacks would result in high signal-to-noise measure- 
ments dominating, but these measurements need not be in 
agreement with each other, because the population has an 
intrinsic dispersion. For example, if one has two quasars with 
fluxes 100 ± lmjy and 0± lmjy, what can one say about the 
average in this population? Clearly, it would not be right to 
quote 50 ± 0.7mJy for this average on the evidence of those 
two quasars. 

We have opted to regard flux measurements of individ- 
ual quasars as attempts to measure the mean of the popula- 
tion, so the dispersion in the population is a noise term on 
these measurements. The noise on any particular estimate 
of the population mean is therefore the quadrature sum of 
the flux error and the population dispersion. 

But what is an appropriate value for this population 
dispersion? We have opted to determine this from our data 
simultaneously with the population mean. If Xi are our mea- 
surements, each with a measurement noise level er;, then our 
data should have the following distribution: 

pl x , o) = 1 exp [ M) l 1 (1) 

where fi is the population mean and oo is the population 
dispersion. We find the maximum-likelihood solutions for /t 
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Figure 5. SWIRE 70 (im and 160 /an fluxes at the positions of SDSS quasars, as a function of the quasar absolute magnitudes and 
redshifts. A background sky level has been subtracted, as discussed in the text. The red line shows the prediction for zero flux, while the 
green curves show the location of an M82 starburst SED normalised to uL v = 10 11 , 10 12 and W 13 Lq at 100 /im. Note that at redshifts 
z < 3, and at most absolute magnitudes, the quasars appear significantly above the zero flux line. Error bars have been suppressed for 
clarity, except for a single representative error bar in each panel. Flux errors were estimated from Gaussian fits to the pixel histograms. 



Wave- 0.5 < z < 1.5 0.5 < z < 1.5 1.5 <z < 2.5 1.5 < z < 2.5 0.5 < z < 2.5 0.5 < z < 2.5 

length no far-IR bright QSOs no far-IR-bright QSOs no far-IR bright QSOs 

70 /mi 7.64±0.95 5.98±0.56 6.11±0.65 5.48±0.56 6.98±0.61 5.77±0.40 

160/mi 19.3±3.0 14.2±2.3 16.1±3.6 12.3±3.0 17.9±2.3 13.4±1.8 

Table 2. Mean far-infrared fluxes in mjy for SDSS quasars in the indicated redshift intervals. The errors quoted are the usual error on 
the mean, a/y/~N, where a is the sample standard deviation and N is the sample size. Quasars are deemed to be far-infrared bright in 
the 70 (im stack if they have a flux of over 20 mjy. In the 160 fim stack, the corresponding threshold is 70 mjy. Sky background levels 
have been subtracted as discussed in the text. 



and o"o from our data by maximising Ilp(xi,ai), and esti- 
mate the parameters' 68% confidence bounds directly from 
the likelihood surface. Note that there is no covariance be- 
tween these parameters; this follows from the fact that the 
expected values of measurements are independent of their 
noise levels. We used numerical simulations to verify that 
our maximum-likelihood solutions using this procedure are 
unbiased estimators of the underlying values, and to verify 
our confidence bounds. 

This estimator worked well where we had more than 
three objects in a bin, but encountered numerical problems 



with some bins containing only two or three objects. There 
is also not enough information to constrain both parameters 
when there is only a single object. One option is to neglect 
the population dispersion ao in these problematic cases, but 
this would lead to over-optimistic error estimates. Another 
option is to ignore these bins altogether. For readers that 
wish to do so, the numbers in each bin will be given. How- 
ever, we found that the quantity in bins with > 3 objects 
ranged from 0.51 to 1.41, with a mean 0.84 and standard 
deviation 0.24. On this basis, we chose to adopt 0.84p as 
our estimator for ctq in bins with three or fewer objects. 
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Figure 6. Histogram of map fluxes (fine bins, black line, left-hand ordinates) compared to the fluxes at the positions of the SDSS quasars 
(broad bins, red line, right-hand ordinates). Note the clear excess of positive fluxes at the positions of the SDSS quasars, i.e. the red lines 
are higher than the black lines on the right hand sides. 



We chose redshift and absolute magnitude bins and 
made a noise-weighted stack of the starburst far-infrared 
luminosities of the quasars in figure [T] in each of these bins, 
using the procedure above. Far-infrared luminosities were 
estimated assuming an M82 SED. The results are listed in 
table |3] and presented graphically in figures [7] and [8] We 
will discuss the effects of varying the SED in section 3.4. We 
used 160 /im fluxes as our far-infrared luminosity estimator 
for the SDSS quasars. For PG quasars we used the 60 /im 
fluxes at redshifts up to 0.3, and 100 (im fluxes at higher red- 
shifts. Our results are not sensitive to the choice of z = 0.3 
for this transition. 

Taking bins in isolation, there are no high signal-to- 
noise detections. Despite this however, there are trends ap- 
parent in figures [7] and [8] Firstly, at all redshifts there ap- 
pears to be a significant optical luminosity dependence, scal- 
ing as L° p 4 t 41±0 - 069 at z < 2 (see figure and table and 
with a shallower scaling at higher redshifts. Secondly, at all 
absolute magnitudes below —28 and redshifts z < 2, there 
is a weaker signature of evolution, (see figure [8] and table [5}. 
The poor \ 2 for the power-law evolution model in the —26 > 
^abs > —28 bin is due to the highest redshift data point, in 
which the positive evolution is reversed, curiously mirroring 
the evolution of quasar number density. If this data point is 
excluded, the power-law fit parameters for this J a b s slice (ta- 
ble EJ are v WO L W o = (0.53±0.25) x 1O 12 L x (l + z) 1A6±0A3 , 
with x 2 = 0.92 and Pr(x 2 , v) = 0.37. This is consistent with 
the evolution at brighter absolute magnitudes, which justi- 
fies an estimate of the average evolution rate of 7 a bs > — 28 
quasars of (1 + ^) 1 ' 57±0 ' 29 . 



3.3 Stacked black hole mass results 

We next examined our data for trends with black hole mass. 
The black hole mass computation is based on the extrapola- 
tion of the reverberation mapping technique, which consid- 
ers the velocity (full-width half maximum) of emission lines 
and relates the size of the Broad Line Region (BLR) to the 



continuum luminosity. Assuming the dynamics of the BLR 
are dominated by the gravity of the black hole, the black 
hole mass is then expressed as 



M BH ^ i?BLR x v 2 /G 



(2) 



where Rblb, is the radius of the BLR and v is the velocity 
of the emission line gas. The velocity v is estimated from 
the FWHM of Hp, Mgll or CIV depending on the redshift 
(Kapsi et al. 2000). For quasars with redshifts greater than 
2 ~ 0.8, is not present in optical spectra, and Mgll and 
CIV have been suggested as alternative estimators. For a 
detailed analysis of the method and the use of the various 
emission lines see Kapsi et al. 2000, McLure & Dunlop 2004, 
and Warner et al. 2004. Depending on the redshift of the 
object and up to a redshift of 4.8, the following relations 
are used, that have been argued by these authors to provide 
equivalent estimates: 



0.0 < z < 0.; 



0.8 < z < 2.1 



2.1 < z < 4.8 



= 4.7 ( 



-V 



M@ 



Mbh 

Ma, 



3.2 



Vio 37 w; 



/ £3000 \ 
uo 37 w; 



= 1a( 



L 



1450 



V 



uo 37 w/ 



/ FWHM(H^) 

V km/s 

/ FWHM(MgH) 

V km/s 



/ FWHM(CIV) V 
I km/s ) 



(3 



where L\ is the luminosity at wavelength A. Black hole 
masses were estimated for objects with available SDSS DR5 
spectroscopy. The FWHM of the lines were derived from the 
emission line sigmas given by the SDSS pipeline, present in 
the second extension of the fits spectra. The error bars on 
the black hole masses are estimated from the uncertainty of 
the FWHM of the lines. We also searched the literature for 
additional black hole mass estimates. For Palomar-Green 
quasars in Boroson & Green (1992), we recomputed the 
black hole masses using the relations above. For the z = 6.28 
quasar SDSS J1148+5251, we applied equation[3]to the data 
in Shields et al. (2006). For other z > 5 quasars, black hole 
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Figure 7. Optical luminosity dependence of the far-infrared luminosities of quasars, in five redshift slices, assuming an M82 SED 
template. Note that the far-infrared luminosities scale with the /-band luminosities, roughly as Lpir <x L®-A . 









0.05 < z < 0.5 


0.5 < z < 1 


1< z < 2 


2 < z < 4 


4 < z < 7 


labs < — 28 








5.5 ±4.5 (2) 


4.9 ± 2.0 (31) 


5.43 ± 0.76 (66) 


4.32 ± 0.70 (63) 


-26 > 7 abs 


> 


-28 


0.52 ±0.38 (9) 


1.9 ± 1.4 (3) 


2.17 ± 0.42 (48) 


3.71 ± 0.87 (76) 


2.41 ± 0.47 (55) 


-24 > / abs 


> 


-26 


0.33 ±0.13 (28) 


0.30 ±0.11 (25) 


0.83 ± 0.22 (74) 


0.7 ± 2.3 (3) 


2.0 ± 2.2 (1) 


-22 > / abs 


> 


-24 


0.086 ± 0.018 (55) 


0.18 ±0.08 (20) 


0.62 ±0.66 (2) 




1.7 ± 1.9 (1) 


4bs > -22 






0.041 ± 0.021 (3) 











Table 3. Error-weighted mean far-infrared vL v luminosities of quasars at 100 /im in redshift and luminosity bins, in units of 10 Lq and 

assuming an M82 SED shape. Numbers in brackets are the number of quasars in each bin. The errors are ( £o\ J , as appropriate 
for error-weighted means, where the errors Oi include an estimate of the population dispersion. 
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Figure 8. Evolution of quasar far-infrared luminosities in four absolute magnitude slices. Note the trends to brighter luminosities at 
higher redshifts, except in the brightest absolute magnitude strip. This can be interpreted as star formation in quasars decreasing with 
cosmic time, except perhaps in the most luminous quasars. 



z range 


Pi 


Vi 




Pr(x 2 , v) 




Pr(x 2 , v) 












(no 7 abs dependence) 


0.05 > z > 0.5 


28.3 ± 1.5 


-0.196 ± 0.057 


0.52 


0.59 


2.67 


0.046 


0.5 > z > 1 


27.00 ± 0.81 


-0.220 ± 0.070 


0.62 


0.53 


1.23 


0.295 


1 > z > 2 


25.21 ± 0.53 


-0.157 ± 0.036 


0.38 


0.68 


3.96 


0.008 


2 > z > 4 


19.4 ± 5.5 


-0.070 ± 0.038 


0.77 


0.38 


2.58 


0.076 


z > 4 


22.1 ± 2.9 


-0.080 ± 0.034 


0.05 


0.95 


1.90 


0.127 



Table 4. Fits to the data in figurc[7] assuming a functional form log 10 (i^iooiloo/lO 12 i0) = P2(^ a bs ~ Pi)- The noise-weighted average 
value of p2 at z < 2 is —0.176 ± 0.028, corresponding to fioo-^ioo °c ^opt 41± ° ° 69 • Also given are the x 2 results for the best model with 
no ^abs dependence. 
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evolution) 


-22 < / abs < 


-24 


0.055 ± 0.015 


1.93 ± 0.57 


0.16 


0.85 


0.89 


0.444 


-24 < / abs < 


-26 


0.180 ± 0.077 


1.43 ± 0.54 


0.60 


0.61 


1.39 


0.234 


-26 < / abs < 


-28 


0.97 ± 0.30 


0.57 ± 0.21 


2.73 


0.04 


4.59 


0.001 


-28 < / abs < 


-32 


7.3 ± 3.5 


-0.26 ± 0.31 


0.23 


0.80 


0.39 


0.759 



Table 5. Fits to the data in figure 151 assuming a functional form fioo£<ioo/lO 12 L0 = P3(l + z) Pi . Also given are the x 2 results for the 
best non-evolving model. The bin —26 < / a bs < —28 has a poor x 2 for the evolving model due to the highest-redshift bin, as discussed 
in section [3.2l Excluding this data point, the average evolution at 7 abs < —28 is (1 + z )!-57±0.29 
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0.05 < z < 0.5 0.5 < z < 1 1 < z < 2 2 < z < 4 4 < z < 7 

10 5 < M BH < 10 6 2.4 ± 3.6 (1) 2.5 ± 2.4 (2) 

10 6 < M BH < 10 7 0.017 ± 0.012 (2) 9.3 ± 8.6 (1) 0.54 ± 0.49 (2) 

10 7 < M BH < 10 8 0.25 ±0.13 (21) 2.4 ±1.8 (3) 2.8 ±1.1 (6) 

10 8 < M BH < 10 9 0.121 ± 0.034 (20) 1.7 ±2.7(1) 5.4 ± 1.0 (11) 2.7 ±1.3 (13) 

10 9 < A/ BH < 10 10 0.37 ±0.35 (3) 10.0 ±9.0(1) 5.0 ±2.7 (14) 2.8 ±1.1 (15) 4.8 ±1.4 (11) 

Table 6. Error-weighted mean far-infrared uL v luminosities of quasars at 100 fim in redshift and black hole mass bins, in units of 

1O 12 L0 and assuming an M82 SED shape. Numbers in brackets are the number of quasars in each bin. The errors are (£cr~ 2 ) , as 
appropriate for error-weighted means, where the errors <Ti include an estimate of the population dispersion. 



masses were taken from Kurk et al. (2007) and Jiang et al. 
(2007), using our equation [4] where possible, and equation[5] 
where only CIV is available. The black hole mass estimates 
vary between these two groups, even when using consistent 
conversions based on the same emission line, but because our 
mass bins are very broad (1 dex), these fractional variations 
(~ 50%) are not enough to move quasars to different bins. 
The results of the stacks in black hole mass and redshift 
bins are given in table [6] and plotted in figures [9] and 1101 
These figures give the host galaxy far-infrared luminosities 
as a function of black hole mass and of redshift. Tables [Jj 
and [8] demonstrate the statistical significances of the trends 
in these figures. 



3.4 Robustness to SED assumptions 

Any multi-wavelength compilation such as this will in- 
evitably rely on SED assumptions to relate the multi- 
wavelength observations. We have assumed an M82 SED up 
to this point, and quoted far- infrared luminosities on that 
basis, but this will inevitably neglect any AGN dust tori 
contributions in the mid-infrared. At best, our far-infrared 
luminosities can only be considered estimates of the star- 
burst bolometric contributions. Starbursts, too, have a vari- 
ety of SEDs, and in this section we will test the robustness 
to our assumed starburst SED shape. The most far-infrared- 
luminous quasars (i.e. those with direct SWIRE detections) 
were found by Hatziminaoglou et al. (2008) to resemble the 
heavily- obscured starburst Arp 220 more often than M82, 
though they conjectured that fainter quasars would be more 
likely to resemble M82. We tested our SED dependence by 
re-running our stacking analyses with an Arp 220 spectrum. 
Reassuringly, very similar trends are present as in the M82 
case, perhaps because most of the rest-frame measurements 
are in spectral regions where the M82 and Arp 220 SEDs are 
similar (see tables [9l and 1 10[) . Moreover, it would not appear 
to be possible to attribute the trends in figures 171151 191 and l 101 
to variations in SED shape as a function of quasar absolute 
magnitude, black hole mass or redshift. 

Few objects in our compilation have photometry at 
more than one wavelength, except the SWIRE quasars. We 
tested our SED assumption by fitting M82 and Arp 220 
SEDs to the SWIRE photometry. Of 236 quasars with pho- 
tometric measurements at both 70 /im and 160 /im, 167 had 
X 2 < 1 for either M82 or Arp 220 SEDs. Of these, 96/167 
(57%) had a lower X 2 for M82. Since most of this photom- 
etry is non-detections, all this is capable of showing is that 
the underlying mean SED is more likely to resemble M82 
than Arp 220, in agreement with the suggestion of Hatz- 
iminaoglou et al. (2008). Objects with \ 2 > 1 typically had 



evidence of a mid-infrared excess, suggestive of an AGN dust 
torus. 

We also tried estimating comparing the 70^im:160^im 
flux ratio with the predictions for redshifted Arp 220 and 
M82 SEDs. The results are shown in figure [TT] If an SED 
model is correct on average, the histogram in figure [TTJ 
should be centred around zero. This again suggests that 
an M82 SED is a better match to the average quasar far- 
infrared-submm SED than Arp 220 (though the most lumi- 
nous may nevertheless more resemble Arp 220). This is in 
contrast with the bright quasars studied by Hatziminaoglou 
et al. (2008), though in keeping with their suggestion that 
fainter quasars have less heavily obscured SEDs. Note how- 
ever that we have already excluded SED dependence on op- 
tical luminosity, redshift or black hole mass as explanations 
for the trends in figures [TJ [8j [9] and [10] 



4 DISCUSSIONS 

4.1 Predictions for Herschel 

Stacking analyses in general only yield information on the 
mean fluxes, and yield little information on the dispersion 
within the population. However, in the case considered here, 
we find evidence for a subset of quasars with bright far- 
infrared fluxes up to ten times the mean fluxes within the 
population. 

The 40 beams/source confusion limits predicted for Her- 
schel by Rowan-Robinson (2001) are 4.6 mjy at 70 /im and 
59m,Iy at 160 /im. If we accept this estimate, then most 
of the targets in the AGN survey would therefore be very 
challenging for PACS direct detection. However, HSPOT re- 
ports a 5cr confusion limit of only 1.2 mjy at 110 /im, so this 
forecast may be pessimistic. 

We used a crude fit to the data in table|3]to predict the 
individual SDSS quasar fluxes (figures [7J and |S), from which 
we estimate that the Herschel ATLAS survey will detect 92% 
of SDSS quasars at z < 0.2 (though the z < 0.2 quasars rep- 
resent only 0.5% of all SDSS quasars). At higher redshifts, 
only the far-infrared-loud subset will be detectable. We esti- 
mate that 66% of SDSS QSOs with far-infrared luminosities 
5x larger than the mean will be detectable in this survey, 
corresponding to about 221(/5/0.05) quasars detected over 
~ 500 deg 2 , where /s is the fraction of quasars with lumi- 
nosities 5x larger than the mean. The detected fraction of 
5x over-luminous quasars at z > 3 is only 5%. However, the 
Herschel ATLAS survey should detect about 98% of all the 
SDSS quasars with luminosities lOx larger than the mean, 
i.e. about 333(/io/0.05) quasars over 500 deg 2 where /io the 
fraction with luminosities 10 x larger than the mean. 

We have neglected type 2 AGN in this analysis. These 
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Figure 9. Host galaxy far-infrared luminosities as a function of black hole mass, for selected bins in table [6] 



2 range 

0.05 > z > 0.5 
2 > z > 4 
2 > 4 



P5 

0.42 ± 0.14 
0.000 ± 0.082 
0.247 ± 0.082 



Pe 

-10.66 ± 0.86 

3 x 10 3 ± 1 x 10 6 

-6.70 ± 0.75 



xl 



1.20 
1.48 
0.87 



Pr(xV) 

0.302 
0.217 
0.454 

12 , 



xl Pr(x 2 .f) 
(no Tabs dependence) 
4.01 0.007 
1.11 0.348 
3.03 0.017 



Table 7. Fits to the data in figure |91 assuming a functional form logjQ^iooiioo/lO 12 !/©) = ps(log 10 (AfBH) +Pe)- Also given are the 
X 2 results for the best model with no Mbh dependence. 



will double or triple the total number of AGN detected by 
the Herschel ATLAS survey. AGN not detected individually 
in this survey will be detectable in stacking analyses. It will 
be illuminating to test whether sub-classes of quasars have a 
greater tendency to be far-infrared-loud in this survey (e.g. 
broad absorption line quasars, nitrogen-rich quasars). 



4.2 Physical interpretation 

It might be possible for AGN dust tori models to account 
for the far-infrared and submm luminosities of quasars, but 
only by assuming very high equatorial optical depths and 
large physical sizes. If quasar heating dominated the far- 
infrared outputs throughout our sample, we would expect 
a linear correlation between quasar absolute magnitude and 
far-infrared luminosity (figure [7}, whereas the observed cor- 
relation is shallower. While we cannot exclude AGN heat- 



ing in a subset of our objects, we will follow Efstathiou 
& Rowan- Robinson (1995) in treating the far- infrared and 
submm luminosities of quasars as being typically dominated 
by star formation. 

The far-infrared luminosities scale linearly with the star 
formation rates as 



SFR 



-Mq /yr 



(6) 



5.8 x 10 9 L Q " 

where Lfir is the bolometric luminosity from the starburst 
(see Kennicutt 1998), assuming a Salpeter initial mass func- 
tion (IMF) from 0.1 to 100 Mg, implying that our quasars 
are forming stars at around 200 — 2000 Mq /yr. The local 
spheroid associated with a 10 8 (10 9 ) Mq black hole has a 
mass of about 5 x 10 10 (5 x 10 11 ) Mq (Marconi & Hunt 2003, 
Hairing & Rix 2004), though there are indications that the 
spheroids are around a factor of 4 less massive at z = 2 (e.g. 
McLure et al. 2006). Sustained star formation at these rates 
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Figure 10. Host galaxy far-infrared luminosities as a function of redshift, for selected bins in table [6] 



M BH range 


P7 


P8 




Pr( X V) 


xl 


Pr(xV) 












(no 


evolution) 


10 6 < M BH < 10 7 


0.0099 ± 0.0085 


2.14 ± 0.71 


1.12 


0.289 


1.14 


0.319 


10 7 < M BH < 10 8 


0.18 ± 0.11 


1.49 ± 0.38 


0.31 


0.581 


3.55 


0.029 


10 8 < M BH < 10 9 


0.082 ± 0.024 


2.12 ± 0.24 


8.47 


0.0002 


10.71 


4 x lO" 7 


10 9 < M BH < 10 10 


0.36 ± 0.26 


1.40 ± 0.43 


1.04 


0.374 


4.09 


0.003 



Table 8. Fits to the data in figure llOl assuming a functional form j/ioo-kioo/lO 12 -^© = P7(l + z) ps . Also given are the x 2 results for the 
best non-evolving model. 



could assemble the z — stellar mass hosts in a few x 10 7 
to a few xlO 9 years. 

The correlation between inferred black hole accretion 
and star formation is similar to the one reported by Hao et 
al. (2007), though they combined low luminosity - low red- 
shift quasars with high luminosity - high redshift quasars, 
so their results could also have been attributable to evolu- 
tion. We span a much bigger range of the optical luminosity 
- redshift parameter space (figure [TJ, so these caveats do not 
apply to our results. 

The quasar magnitudes in our lowest redshift bin may 
in principle be contaminated by the host galaxies. If so, cor- 
recting for this effect would only strengthen the dependence 
of star formation on quasar luminosity, since the correc- 
tion would apply preferentially at the faintest optical mag- 
nitudes. 



One difficulty in the interpretation of these results is the 
possibility of luminosity-dependent reddening of quasars. 
However, even an Ay of one at the faintest end and zero at 
the brightest would have little impact on our correlations, 
given the size of our errors and range of absolute magni- 
tudes. We have assumed a single optical spectral index for 
our quasars, but a slightly better approach would be to use 
the optical spectra themselves, correcting for dust absorp- 
tion using the Balmer decrement where possible, or using 
rest-frame hard X-ray luminosities. Alexander et al. (2005) 
found an approximately linear relationship between hard X- 
ray and far-infrared luminosities in a heterogeneous sample 
of AGN-dominated submm-selected galaxies, though if one 
adds the submm galaxies classified as starbursts, their cor- 
relation is shallower with a wide dispersion. We have not 
excluded candidate starburst-dominated objects. Alexander 
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0.05 < z < 0.5 0.5 < z < 1 1 < z < 2 2 < z < 4 4 < z < 7 

4bs < -28 8.4 ±6.9 (2) 6.2 ± 2.8 (31) 3.95 ± 0.52 (66) 3.96 ± 0.64 (63) 

-26 > I &hs > -28 0.53 ±0.38 (9) 2.3 ±1.6 (3) 2.51 ± 0.48 (48) 4.4 ± 1.1 (76) 2.19 ± 0.43 (55) 

-24 > / abs > -26 0.47 ± 0.20 (28) 0.29 ± 0.11 (25) 0.84 ± 0.21 (74) 1.0 ±3.5 (3) 2.0 ±2.2(1) 

-22 > 7 abs > -24 0.107 ±0.021 (55) 0.18 ± 0.08 (20) 0.59 ± 0.63 (2) 1.7 ±2.0(1) 
Jabs > -22 0.055 ± 0.028 (3) 

Table 9. Error-weighted mean far-infrared vL v luminosities of quasars at 100 /im in redshift and luminosity bins, in units of 10 12 Z/q and 

assuming an Arp220 SED shape. Numbers in brackets are the number of quasars in each bin. The errors are (Sov j , as appropriate 
for error-weighted means, where the errors Oi include an estimate of the population dispersion. 



0.05 < z < 0.5 0.5 < z < 1 l<z<2 2 < z < 4 4 < z < 7 

10 5 < M BH < 10 6 1.8 ±2.7(1) 2.4 ±2.3 (2) 

10 6 < M BH < 10 7 0.023 ± 0.017 (2) 7.2 ±6.6(1) 0.46 ± 0.42 (2) 

10 7 < M BH < 10 s 0.35 ±0.19 (21) 1.6 ±1.1 (3) 2.38 ± 0.90 (6) 

10 8 < M BH < 10 9 0.183 ± 0.053 (20) 1.0 ±1.6(1) 4.25 ± 0.80 (11) 2.5±1.2(13) 

10 9 < M BH < 10 10 0.37 ±0.34 (3) 15 ± 14 (1) 6.0 ±3.8 (14) 2.16 ± 0.86 (15) 4.6 ±1.3 (11) 

Table 10. Error- weighted mean far-infrared vL v luminosities of quasars at 100 ^m in redshift and black hole mass bins, in units of 

10 12 L© and assuming an Arp220 SED shape. Numbers in brackets are the number of quasars in each bin. The errors are (Scr. ) , 
as appropriate for error-weighted means, where the errors it; include an estimate of the population dispersion. 



et al. (2005) also demonstrate that local AGN show a large 
scatter in their star formation - black hole accretion relation- 
ship; our lowest-redshift quasars in figure [7] show marginal 
evidence for a steeper correlation than at higher redshifts, 
broadly in agreement with these observations of local active 
galaxies. 

The implicit correlations in figure [7] between star forma- 
tion rate and black hole accretion rate hint at common phys- 
ical parameters (such as gas supply) , despite the disparity of 
spatial scales, but in keeping with qualitative expectations 
from the black hole - spheroid connections in local galax- 
ies (e.g. Maggorian et al. 1998, Ferrarese & Merritt 2000). 
There is no insight to be gained by supposing the only com- 
mon parameter is the total mass of the system, i.e. that this 
is simply reflecting only a size dependence, because one must 
still hypothesize some mechanisms to tie both parameters to 
the total mass (e.g. Serjeant et al. 1998); in any case, the 
trends in figure [7] follow approximately I/pir oc L°p t rather 
than a linear relationship. Furthermore, the lack of any obvi- 
ous correlation with black hole mass at redshifts 0.5 < 2 < 4 
(see below) argues against any simple scaling with the size 
of the system, at least outside the local Universe. 

This non-linear relationship and its evolution do not fol- 
low expectations from some semi-analytic models. Accord- 
ing to the model of Croton (2006), the ratio of black hole 
accretion rate and star formation rate is constant with scale 
but increases with redshift. This is partly due to increased 
disk disruption at high redshifts generating starbursts but 
not black hole accretion, and partly due to evolution in the 
black hole feeding rate in this model. The model does suc- 
ceed in reproducing the evidence for evolution in the black 
hole - bulge mass relationships and (qualitatively at least) 
our evolving normalisation of the black hole accretion - star 
formation relation from z < 0.5 to 1 < 2 < 2. However, 
it does not reproduce our scale dependence. Our observa- 
tions may prove to be a useful constraint on AGN feedback 
models. 

If AGN feedback directly regulates stellar mass assem- 
bly in the host, then we may expect stronger trends of far- 
infrared luminosity with black hole accretion than with black 



hole mass. Not all of our sample have black hole mass esti- 
mates, so our tests of dependence on black hole mass and 
redshift are more noisy than our correlations against quasar 
luminosity, though we span a larger logarithmic range of 
black hole mass than quasar luminosity. In the local Uni- 
verse, we see a hint of a relation between star formation 
rate and black hole mass (figure [9] table [Jj. At these lower 
redshifts, the hosts may already have assembled a large frac- 
tion of their 2 ~ stellar masses, so this may partly repre- 
sent mutual size dependence. However, at 0.5 < 2 < 4 there 
seems to be no evidence for a dependence on black hole mass 
(figure [5]) despite the dependence on inferred black hole ac- 
cretion (figure[7J. At 2 > 4 there is some evidence for a weak 
trend with black hole mass (see table |HJ, varying roughly as 
SFRoc M B (j 4 ; the weakness of this trend suggests that it is 
less closely related to the primary underlying physical mech- 
anism than the SFR-L opt relation. There are nevertheless 
hints of trends with redshift at fixed black hole masses (fig- 
ure |S] tables |S] and [SJ). More black hole mass estimates are 
needed to improve the statistics, but our measurements are 
consistent with feedback from black hole accretion at 2 > 1 
regulating the stellar mass assembly in their hosts. 

It is likely that the e-folding timescale for black hole 
growth (tbh — 4A(0.1/e) x 10 7 yr, where A is the Eddington 
ratio and e the accretion efficiency) is much faster than the 
stellar mass assembly timescale (e.g. Malbon et al. 2007). 
Typical quasar lifetimes at 2 < 3.5 may not be much longer 
than a single e-folding scale (e.g. Martini & Weinberg 2001, 
Shen et al. 2007) , though may be several e-foldings at higher 
redshifts (Shen et al. 2007). There are 3.6 e-foldings from 
7ab = —22 to Jab = —26, making it unlikely that the re- 
lationships in figure [7] represent the evolutionary tracks of 
individual objects. 

We have shown that active nuclei are on average as- 
sociated with luminous or ultraluminous starbursts at all 
redshifts, and all absolute magnitudes brighter than about 
Tab = —22. This relationship does not on its own help us ad- 
dress whether the AGN initiates the starburst or is started 
concurrently, or whether the AGN occurs at some midway 
point, or whether the AGN quenches the starburst. How- 
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Figure 11. Comparison of the 70/jm:160 /an flux ratio for SWIRE SDSS quasars (including individual SWIRE non-detections) for an 
M82 SED (thick histogram) and an Arp 220 SED (thin histogram). The left panel is restricted to objects at redshifts z < 2, and the right 
panel has the additional restriction that the best \ 2 °f the M82 and Arp 220 models is better than 1. The x-axis shows the observed 
flux ratio divided by the predicted ratio, minus one. If the SED model choice is correct on average, the histogram should be centred on 
zero. The M82 histogram (thick) appears better centred on zero, so is favoured by this test. 



ever, if quasar lifetimes are as long as 600 Myr at z > 3.5, 
which is the upper limit to the lifetimes suggested by Shen 
et al. (2007), then we would expect AGN feedback to have 
quenched the star formation in nearly all z > 3.5 quasars. 
Our tentative high-z detections suggest this is not the case. 
The high-z constraints will shortly be made much stronger 
with the large far-infrared and submm photometric surveys 
of z > 3.5 quasars from the Herschel ATLAS key project, 
SCUBA-2 and other facilities. Conversely, the shorter in- 
ferred quasar lifetimes at lower redshifts, the lack of evidence 
for any dependence of star formation on black hole mass, the 
observed dependence of star formation rate with quasar lu- 
minosity, and the local bulge - black hole relationships, are 
all consistent with feedback from black hole accretion regu- 
lating stellar mass assembly at lower redshifts. 
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APPENDIX A: DISCUSSION OF 
PALOMAR-GREEN QUASARS 

We use our SCANPI 100 fim measurements in preference 
to other IRAS 100 /im measurements (e.g. Saunders et al. 
1989), but that withstanding, we use the 100 /im photome- 
try with the smallest errors with the exceptions of the cases 
discussed below. At 60 fim we use the measurements with 
the smallest errors regardless of their source, again with the 
exceptions discussed below. Haas et al. (2000, 2003) do not 
quote errors on their ISO photometry, but state that the 
detections range from 3— 10cr. We conservatively assume 3a 
unless stated otherwise. Objects dominated by non-thermal 
emission at 60-100 fim have been eliminated from our stack- 
ing analyses. The adopted photometry for the Palomar- 
Green sample is given in table IAT1 

0007+106, or MRK1501, is a radio-loud flat-spectrum 
quasar, and is likely to have non-thermal emission dominat- 
ing at 60-100 fim. 

0050+124, or lZwl, has a 60 (100) /im measurement 
from ISO of 1752 (2339) mjy reported by Haas et al. 2000, 
but these are inconsistent with our SCANPI measurements 
of 2161+52 (1749+187) mjy which show no obvious anoma- 
lies in the fits. We have opted to use our SCANPI photom- 
etry. 

0157+001 has a 60 /im flux measurement of 2210 mjy 
reported by Haas et al. 2003, though errors are not 
quoted. Sanders et al. 1989 quotes an IRAS measure- 
ment of 2377+56 mjy. Our IRAS SCANPI measurement of 
2348+73 mjy is consistent with the latter rather than the 



former. We have opted to use our SCANPI photometry, 
which shows no obvious anomalies in the fit. 

0832+251 (z=0.320) is reported as < 126 mjy at 60 ^m 
in Sanders et al. 1989 but has a SCANPI measurement of 
352+60 mjy. There is a similar discrepancy at 100 /im. How- 
ever this is because SCANPI's fit is strongly affected by the 
nearby IRAS galaxy IRAS 08325+2512 at z=0.017, which is 
2.5 arcmin from the QSO. Both the 60 and 100 /im coadded 
scans appear to be fairly flat off-source. We therefore found 
the maximum-likelihood fit for the amplitudes of a source 
fixed at the target position, and another with a position 
allowed to vary. 

0838+770 has an ISO 100 /xm flux measurement of 
180 mjy from Haas et al. 2003, which disagrees with the 
Sanders et al. 1989 IRAS measurement of 426+30 mjy. Our 
IRAS SCANPI measurement is 293+184 mjy, though with 
strong baseline drifts in the coadded timeline. Given the un- 
certainties we have uncovered in the IRAS 100 /im flux cali- 
bration at this level, and the baseline drifts in our SCANPI 
data, we have opted to use the Haas et al. photometry with 
an assumed 33% error. 

1001+054 has an ISO 60 /im flux measurement of 
140 mjy from Haas et al. 2003, which disagrees with our 
IRAS SCANPI measurement of 23+49 mjy. There are no 
hints of flux at the target position in our coadded scans. 
Sanders et al. 1989 report 27+9 mjy, which has a remark- 
ably small quoted error. Nothing is reported at this position 
in either the IRAS PSC or FSC, and nothing is evident on 
the ISSA plates. We have opted to use the Sanders photom- 
etry. 

1022+519 has a 100 /im flux measurement from the 
IRAS Faint Source Reject catalogue of 798+176 mjy. Our 
SCANPI photometry is 200+103 mjy, with a fairly stable 
coadded baseline. Owing to this, and the lower quoted error 
of our SCANPI measurements, together with the uncertain- 
ties we have uncovered in the catalogued IRAS 100 micron 
fluxes at this level, we have opted to use our SCANPI mea- 
surement. 

1100+772, or 3C249.1, is radio-loud and probably 
synchrotron-dominated at both 60 fj,m and 100 /im. 

1103-006 has a 60 [im IRAS flux of 130+51 mjy quoted 
in Sanders et al. 1989, while our SCANPI measurement is - 
8+98 mjy, though the coadded scans are affected by baseline 
drifts. We have used the Sanders measurement. 

1211+143 The Haas et al. 60 /im measurement of 
518 mjy disagrees both with the Sanders et al. 1989 mea- 
surement of 305+53 mjy and our own IRAS SCANPI mea- 
surement of 284+81 mjy. At 100 /im the disagreement is 
more striking, with Haas et al. reporting an upper limit of 
< 279 mjy, while our SCANPI measurement is 427+182 mjy 
and Sanders et al. report 689+119 mjy. Although baseline 
drifts are evident in our coadded scans, the background sub- 
tractions at the position of our target appear to be reliable, 
and the profile is well-fit. We have opted to use the lowest 
noise IRAS measurements. 

1226+023, or 3C273, has a clearly non-thermal spec- 
trum in IRAS passbands. 

1302-102 is radio- loud and probably synchrotron- 
dominated at 60 and 100 /im. 

1351+640 has a 100 /an flux measurement from Haas et 
al. 2003 of 526 mjy, but the IRAS SCANPI measurement is 
912+156 mjy. Although baseline drifts are apparent in our 
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coadded scans, the background subtractions at the position 
of our target appear to be reliable. Sanders et al. 1989 report 
1184±26mJy. We have opted to use our SCANPI measure- 
ment. 

1501+106 has a 60 /jm measurement from Haas et 
al. 2003 of 750 mjy, but our SCANPI measurement is 
473±36mJy. Sanders et al. 1989 report 486±42 mjy. We 
adopt our SCANPI measurement as the most likely lowest- 
noise choice. 

1545+210 is radio-loud and probably synchrotron- 
dominated at 60 and 100 fim. 

1613+658 has a 100 /jm measurement from Haas et al. 
2000 of 1002 mjy, consistent with the Sanders et al. 1989 
measurement of 1090+59 mjy, but our SCANPI measure- 
ment is 474+70 mjy. Our SCANPI measurement shows a 
baseline drift that may be over-corrected, so we opt to use 
the Haas et al. measurement, and assume an error of 10%. 

1704+608, or 3C351, is radio-loud and probably 
synchrotron-dominated at 60 and 100 fim. 

1718+481 is radio-loud and probably synchrotron- 
dominated at 60 and 100 /im. 

2209+184 is radio-loud and probably synchrotron- 
dominated at 60 and 100 jj,m. 

2251+113 is radio-loud and probably synchrotron- 
dominated at 60 and 100 /im. 

2344+092 is radio-loud and probably synchrotron- 
dominated at 60 and 100 /im. 

2349-014 is radio-loud and probably synchrotron- 
dominated at 60 and 100 fim. 

Only contains objects which are not dominated by syn- 
chrotron at 60-100um 
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00 
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00 


53 
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54 
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01 
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0. 
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01 


59 


50 
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23 40.62 


2348+73 
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0. 
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08 


10 
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1 


0832+251 


08 
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35 
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0838+770 
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.26 
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0. 


029 
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11 


21 


47 
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+11 44 18.26 


452+49 


481+267 


0. 


049 


1121+422 


11 


24 


39 
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11 


29 
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16 
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Table Al. This table lists the adopted photometry for the Palomar-Green quasar sample. The columns give the name, the position, the 
60 jim and 100 //m photometry, and the redshift. (Continued overleaf.) 
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Name Right Ascension Declination Sqq SlOO Redshift 

(J2000) (J2000) (mjy) (mjy) 
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652+21 


793+87 


0. 


.079 


1444+407 


14 


46 


45 


.940 


+40 


35 


05. 


70 


57+30 


80+27 


0. 


.267 


1448+273 


14 


51 


08 


.8 


+27 


09 


27 




117+37 


-34+100 





.065 


1501+106 


15 


04 


01 


.201 


+10 


26 


16. 


15 


473+36 


77+144 





.036 


1512+370 


15 


11 


13 


.042 


+36 


50 


50. 


11 


61+20 


160+159 


0. 


.37 


1519+226 


15 


21 


14 


.2 


+22 


27 


43 




-21+49 


155+150 


0. 


.137 


1522+101 


15 


24 


24 


.6 


+09 


58 


30 




37+71 


-38+201 


1. 


.321 


1534+580 


15 


35 


52 


.361 


+57 


54 


09. 


21 


140+51 


136+128 





.03 


1535+547 


15 


36 


38 


.361 


+54 


33 


33. 


21 


61+32 


81+148 


0. 


.038 


1538+477 


15 


39 


34 


.8 


+47 


35 


31 




97+39 


107+121 


0. 


.770 


1543+489 


15 


45 


30 


.240 


+48 


46 


09. 


10 


348+26 


371+79 





.4 


1552+085 


15 


54 


44 


.6 


+08 


22 


22 




-82+103 


-276+140 


0. 


.119 


1612+261 


16 


11 


13 


.210 


+26 


04 


16. 


20 


252+72 


330+629 


0. 


.131 


1613+658 


16 


13 


57 


.179 


+65 


43 


09. 


.58 


635+19 


1002+100 


0. 


.129 


1617+175 


16 


20 


11 


.288 


+17 


24 


27. 


70 


52+45 


45+108 


0. 


.114 


1626+554 


16 


27 


56 


.0 


+55 


22 


31 




-28+46 


70+23 





.133 


1630+377 


16 


32 


01 


.120 


+37 


37 


50. 


00 


5.9+36 


-105+110 


1. 


.466 


1634+706 


16 


34 


28 


.884 


+70 


31 


33. 


04 


318+23 


444+80 


1. 


.334 


1700+518 


17 


01 


24 


.800 


+51 


49 


20. 


00 


480+36 


374+125 





.292 


1715+535 


17 


16 


35 


.5 


+53 


28 


15 




3.2+60 


-88+123 


1. 


.920 


2112+059 


21 


14 


52 


.6 


+06 


07 


42 




105+19 


193+370 


0. 


.466 


2130+099 


21 


32 


27 


.813 


+10 


08 


19. 


16 


479+12 


485+162 


0. 


.062 


2214+139 


22 


17 


12 


.26 


+14 


14 


20. 


9 


337+11 




0. 


.067 


2233+134 


22 


36 


07 


.680 


+13 


43 


55. 


30 


80+68 


-647+302 





.325 


2302+029 


23 


04 


45 


.0 


+03 


11 


46 




130+66 


118+174 


1. 


.044 


2304+042 


23 


07 


02 


.9 


+04 


32 


57 




60+63 


70+130 


0. 


.042 


2308+098 


23 


11 


17 


.758 


+10 


08 


15. 


46 


83+87 


-539+420 





.433 



Table A2. (continuation of table lAlh This table lists the adopted photometry for the Palomar-Green quasar sample. The columns give 
the name, the position, the 60 (im and 100 (im photometry, and the redshift. 



